N84-22566 


NASA Technical Memorandum 83639 


Determination of Compressor In-Stall 
Characteristics from Engine 
Surge Transients 


Carl F. Lorenzo, Francis P. Chiaramonte, 
and Charles M. Mehalic 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

Twentieth Joint Propulsion Conference 
cosponsored by the AIAA, SAE, and ASME 
Cincinnati, Ohio, June 11-13, 1984 




E-2082 


DETERMINATION OF COMPRESSOR IN-STALL CHARACTERISTICS 
FROM ENGINE SURGE TRANSIENTS 


Carl F. Lorenzo, Francis P. Chiaramonte, 
and Charles M. Mehalic 
National Aeronautics Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Abstract 

A technique for extracting the in-stall 
pumping characteristics for an axial flow 
compressor operating in an engine system 
environment is developed. 

The technique utilizes a Hybrid computer sim- 
ulation of the compressor momentum equation into 
which actual transient data are used to provide 
all terms but the desired compressor characteris- 
tic. The compressor force characteristic as a 
function of corrected flow and speed result from 
the computation. The critical problem of data 
filtering is addressed. 

Results for a compressor operating in a 
turbofan engine are presented and comparison is 
made with the conventional compressor map. The 
relationship of the compressor surge characteris- 
tic with its rotating stall characteristic is 
explored. Initial interpretation of the measured 
results is presented. 

Nomenclature 

A flow area 

F corner frequency, force 

G ( s ) filter transfer function 
g gravitational constant 

K constant 

£ compressor length 

N rotor speed 

P pressure 

PR pressure ratio, P 3 T/P 2 T 

s laplace parameter 

sgn sign of 

T temperature 

t time 

V velocity 

'V volume 

Cl flow rate 

0 density 

t time constant 


Subscripts 




c 

compressor 

i 

input 


CHAR 

characteristic 

m 

map 


cs 

control surface 

M 

normal 


c v 

control volume 

nom 

nominal 


corr 

corrected 

0 

output 


d 

differentiator 

S 

static 


D 

dynamic, derivative 

T 

total 


ex 

external 

w 

wall 


f 

friction 

X 

axial 


f 

friction linear 

2 

compressor 

inlet 

fs 

friction square law 

3 

compressor 

exit 


Introduction 

The understanding of turbomachinery com- 
pressor dynamic behavior is central to the stabil- 
ity modelling of turbine engine systems. This is 
especially true when considering the stagnation 
stall or hung stall phenomena manifested in some 
current engines.^'' In this situation a stall 
occurring in the compressor system is difficult or 
impossible to clear with normal control actions. 

The ability to predict the stall stagnation 
or recovery processes requires detailed under- 
standing of compressor dynamic behavior both in 
the stable (unstalled) zone and considerably off- 
design, that is, in-stall and even into the com- 
pressor negative flow regime. 

Numerous investigators have developed dynamic 
modelling procedures for axial flow compres- 
sors. I 2- '*! Typically these models are intended 
for use in full system models and are premised on 
actuator disk concepts. These concepts combine 
compressor steady state performance maps with the 
conservation equations written in a one dimen- 
sional form. 

Usually in this type of modelling the momen- 
tum equation is written across either a stage or 
the complete compressor. The terms which go into 
the compressor momentum equation are the 
pressure-area terms, the flux of momentum terms, 
and the external force which is generated by the 
pumping action of the compressor. This is the sum 
of the internal forces applied to the compressed 
fluid. It is this forcing function supplied by 
the compressor in the dynamic sense which is to be 
determined. In the past, in order to make the 
momentum equation consistent with steady state 
data, the normal compressor map, that is pressure 
ratio versus flow (or a phi-psi map), was used to 
derive the force applied by the compressor to the 
fluid. Determination of this force in the dynamic 
sense is quite difficult. 

The subject of this paper is to attempt to 
imply this instantaneous force based on experimen- 
tal data. It is the premise of this paper that if 
all the terms required for the momentum equation 
with the exception of the instantaneous force 
applied to the fluid by the blades can be mea- 
sured, then the instantaneous force can be implied 
by the momentum equation. This in turn can be 
related back to the desired compressor character- 
istic. A concept similar to this has beep for- 
warded in an unpublished paper by Gamache'^/ in 
which some initial attempts to determine the com- 
pressor characteristics were implemented. This 
effort resulted in double valued compressor char- 
acteristics. 



The present paper attempts to establish a 
computational technique for determining the com- 
pressor instantaneous force from experimental 
data, to relate that force to the compressor 
in-stall characteristic, and to establish guide- 
lines for data filtering. The paper will present 
typical results for data generated on a turbofan 
engine, examine the effects of compressor speed, 
examine the effects of filtering, and study the 
transition into rotating stall. 

Analytical Basis 


= P 2T A 2 - P 3T A 3 + F X + 7 p 2 A 2 V 2 

- * p 3 A 3 v2 (6) 

Neglecting the residual 1/2 pAV^ terms at low 
Mach number gives 

dW 

= P 2T A 2 " P 3T A 3 + F x (7) 


The analysis assumes an overall compressor 
model, that is, the conservation equations can be 
written across the entire compressor. It is, of 
course, also possible to consider the compressor 
on a stage-by-stage basis assuming that measure- 
ments can be made for all the inputs required as 
described above. 


as the working form of the momentum equation. 

Note it is possible to relax some of the above 
assumptions if both static and total pressures are 
used, this is derived in Appendix A. Solving for 
the external instantaneous force acting on the 
fluid gives 


The one dimensional momentum equation for a 
fixed control volume can be written in the fol- 
lowing form (for example, see ref. 6). 

£ F ex ■ / it * / "Vx "out 

cv cs 

- / pV„V dA. (1) 

-/ N x in ' 

cs 

The external forces in this equation are as 
shown in Fig. 1, include the pressure-area terms 
acting on the front and rear faces of the compres- 
sor fluid, and a net force (called the instantane- 
ous force) F x is applied to the fluid by the 
blades, wall, etc. For the compressor configura- 
tion, assuming one dimensional planar flow, i.e., 
V N = v x and noting 



cv 


and 

/•Vx « - < 3 > 

cs 

Then the momentum equation becomes 

. dW 

Z F ex ■ ~ P 3S«3 * F x -|dT 

+ °3 A 3 V 3 ” °2 A 2 V 2 ^ 

Now for low Mach numbers, total pressure can be 
approximated as 


F x ~ P 3T A 3 ^T^ + g dt 

Some investigators notably Grietzer,( p ) 
hypothesize that at least in the transition from 
surge to rotating stall, additional dynamics occur 
between the steady-state characteristic and the 
instantaneous force. For the rotating stall case 
this is the time required for the stall cell 
growth to occur. A simple first order lag is 
assumed to represent these dynamics, thus: 


• dt 


+ F = F 
x m 


(9) 


determines the instantaneous force from the so 
called steady-state force F m . A block diagram 
for the compressor model with this dynamic element 
assumed to be a first order lag is shown in 
Fig. 2. It is not clear that t c has a value 
other than zero for a pure surge condition. 

(Surge data is to be analyzed.) 

Assuming the above model applies, the instan- 
taneous force F x can be implied from Eq. (8). 

Two assumptions are implicit in this calculation, 
these are: (a) the flow is predominately one- 

dimensional (uniform and planar) i.e., without 
rotating stall or secondary flows, and (b) that 
the flow throughout the compressor can be -approxi- 
mated by a single measurement, i.e. integration 
over the control volume is not necessary. If 

P 3T 5 P 2T’ and W c can be actually measured 
as functions of time for a compressor operating in 
surge, then F x can be determined as shown in 
the simplified block diagram of Fig. 3. Also, in 
this figure, the inverse dynamics of Eq. (9) is 
applied to F x to determine F m . Both F x 

and F m are shown plotted against W c cor „ 
which is assumed to be the basis of correlation. 


P T 



Thus the momentum equation becomes 


The actual calculation to determine the 
forces F x and F m is more involved than 
Fig. 3 indicates because: (a) it is possible that 
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the real data is not pure surge but is contami- 
nated by either rotating stall or other multi- 
dimensional flow effects, (b) the real measure- 
ments can be contaminated by instrument or 
recording noise, (c) pure differentiations of the 
data tend to amplify the high frequency content 
(noise) and must, therefore, be approximated (fil- 
tered) and finally (d) several measurements may be 
required to replicate each of the signals that 
serve as inputs to the calculation. 

With respect to the rotating stall issue, if 
rotating stall is present during the surge its 
effect may sometimes be removed from the data by 
careful filtering. This requires low-pass filters 
to be applied to the signals representing the mea- 
sured variables prior to manipulation in calculat- 
ing F x and F m . These filters, if used, 
must have sufficient bandwidth to allow the 
highest harmonics present from the surge phenome- 
non to pass without significant loss. Simultane- 
ously the filter corner frequency must be suffi- 
ciently below the rotating stall frequency to 
provide attenuation. Filters are also required to 
eliminate or minimize instrument noise and as a 
part of the differenti ation process. The subject 
of filter placement and selection will be dis- 
cussed in detail in a later section. 


determined by examining Eqs. ( 8 ) and (9). At 
steady-state conditions the derivative terms in 
these equations become zero, thus, F m = F x 
and. 


F x “ P 3T A 3 P 2T A 2 " P 2T A 2 



( 11 ) 


Hence the normal compressor map characteristic (at 
constant speed) can be expressed as 



( 12 ) 


Here the subscript CHAR is added to differen- 
tiate between the characteristic and the dynamic 
pressure ratio. In actual practice, only very 
small differences in (P 3 T/P 2 T)CHAR are seen 
between taking PpjAp at its steady value or 
dynamic value in the division P X /P?T^ 2 ‘ 

• 

It can be seen from these equations that, the 
F x characteristic is only a scale factor and 
constant translation different from the compressor 
map characteristic. 


Relative to item D) above for a one- 
dimensional reversing flow field two measurements 
are required for each total pressure; i.e., for- 
ward and aft facing. Thus as flow direction 
switches, the measurement read as total pressure 
must also switch. Also, corrected flow is com- 
puted from a Tee probe or Mach probe using appro- 
priate flow calibrations. 


Combining the above considerations, the com- 
putation for the instantaneous force, F x , and 
the map force, F m , for the most general case 
will appear as shown in Fig. 4. In this block 
diagram, the input filters, G-j(s), following the 
measurement have the function to minimize the con- 
tamination by noise and effects of rotating 
stall. The output filters, G 0 (s), are always 
chosen as matched pairs and serve to help elimi- 
nate noise on the plotted signals. Also, practi- 
cal differentiators have associated with them, 
built-in filters to limit amplification of high 
frequency noise. All filters used in this study 
are first order lags, i.e.. 


G(s) 


Xout(s) 1 

Xin(s) ” ts + 1 


( 10 ) 


where s is the Laplace parameter, and t is 
the filter time constant. The values chosen for 
the filter time constants will be discussed in a 
separate section. 


It is the basic block diagram. Fig. 4, which 
is implemented on a hybrid computer and used to 
determine the compressor characteristic. The 
required input dynamic signals are recorded on FM 
magnetic tape, and trunked directly to the 
computer. 


Steady-State Rel ationships 


The relationship between the force character- 
istic F x or F m , and the normal compressor 
map, pressure ratio versus corrected flow, may be 


Nature of Instantaneous Force 


The instantaneous force F x (fig. 5) can 
be thought of as being composed of three parts 
which are: (1) F w , the force of the wall on the 

fluid, (2) Ff, the frictional force acting on 
the fluid and (3) Fg the dynamic force applied 
to the fluid by virtue of the blade rotation. 

These forces may be considered as analogous to 
electrical voltages in a circuit, that is; Fg 
is similar to a generated force, Ff is analo- 
gous to a resistive force and F w , may be 
thought of as a back EMF or reactive force. Thus, 


F 

x 


= f d- 


(13) 


Now the wall force may be approximated by 


2 ( P 2S + P 3S^ ( A 2 " A 3^ 


For 1 inear friction 


(14) 


fi 


K, W 
fa. c 


(15) 


or for square law friction 

F = K, sgn (W ) W 2 (16) 

fs fs 3 ' c c 

For the linear friction case and using Eq. 8 , the 
blade force, Fg may be found as 


F D - P 3T A 3 P 2T A 2 + 2 ( P 2S + P 3S^ ^ A 2 V 


o dW r 

+ K, W + - -rr^- (17) 

fi c g dt v 


Or, if the differences between total and static 
pressure is ignored and letting A = (A 2 + A 3 ) / 2 
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F - P I 

D ~ 3T 


dW 

P 2T S + K f* * c + J-af 


-3F ( 18 ) 


If Kf^ is known or can be estimated 
these equations (eq. (17) or (18)) can be used to 
gain an insight to the nature of the instantaneous 
dynamic force. This has not been done in the cur- 
rent study. 


By far most of the cycles seen in this transient 
are of the mode I type. Mode II, surge instabil- 
ity, is seen at the beginning of the transient, at 
the transition from surge into rotating stall and 
only very sporadically during the long transient. 
The frequencies of both modes increase by about 
0.5 Hz over the duration of the run. 

Instrumentation and Data Acquisition 


Experimental Aspects 

The data used in this study result from a 
single test run performed on a multi-stage axial 
flow compressor in a turbofan engine which was 
tested in the Propulsion System Laboratory of the 
NASA-Lewis Research Center. Under altitude condi- 
tions (inlet pressure of approximately 1/2 atmo- 
sphere) the normal control for the engine was 
over-ridden and a large triangular fuel flow pulse 
was introduced into the combustor. 

This pulse forced the engine into a surge- 
like instability which was allowed to continue 
(approx 32.5 sec) until the compressor settled 
into rotating stall (a stagnation stall) condi- 
tion. An overview of this transient is shown in 
Fig. 6 where compressor discharge static pressure 
P 35 , interturbine temperature ITT, core mechani- 
cal speed N c and power lever angle PLA are 
plotted versus time. The fuel pulse which causes 
the surge is not shown but lasts about 1/2 second 
into the transient after which normal control is 
resumed. 

That the compressor is in rotating stall 
after t = 32.5 seconds is determined by power 
spectral analysis of data from circumferentially 
located pressure transducers at a common axial 
location (compressor inlet) which indicate appro- 
priate phase differences. The rotating stall 
frequency is about 0.41 times the rotor speed, 
starting at about 91.3 Hz and decreasing with time 
and speed. Steady state and dynamic pressure and 
temperature instrumentation were available both 
upstream and downstream of the compressor. 

Interstage wall static pressures were also 
measured. Some of the internal pressures in the 
compressor during the surge process are shown in 
Fig. 7. Here the transition from stable flow into 
surge (at approx 370 msec) is presented. From 
this figure it appears that the stall originates 
near the 8 th stage based on the pressure increase 
in all upstream stages at stall inception. Also 
the character (and frequency) of the first three 
stall cycles are different from the following two 
cycles. The frequency associated with the first 
stalls is about 10.7 Hz while that for the next 
two (and most of those that follow) is about 5.2 
Hz. During the time periods A and B of Fig. 7, 
all the engine measurements are steady (not time 
varying) except those in the fuel control. It is 
assumed from this that the engine is essentially 
"recovered" but that a control instability (called 
a mode I instability in the remainder of the 
paper) forces a repeated stal 1 /recovery cycle. 

The first cycles (10.7 Hz) are considered to be a 
classical surge, and will be called a mode II 
instability. 


In order to achieve high dynamic response 
pressure measurements, miniature pressure trans- 
ducers were used. These solid state transducers 
are subject to significant zero shift with temper- 
ature. Thus, the instruments were water cooled 
and separate, stable lower response instruments 
are used to measure the steady (DC) level. The 
pressure signals analyzed are composites which are 
the sum of the DC signal, taken at a steady period 
just prior to the transient and the high response 
signal less its steady part at the same reference 
time. In spite of these precautions small addi- 
tional DC offsets are found to exist in the data. 
To account for this, at a single speed (N = 11%) 
minor corrections are made to the DC level of the 
signals so that the recovered periods A and B of 
Fig. 7 agree with a corresponding steady period on 
the normal compressor map. High response, 0.003 
inch diameter thermocouples were used for tempera- 
ture measurements. Corrected flow was measured 
with specially calibrated fore and aft facing 
total pressure (TEE) probes. 

All dynamic measurements were recorded on 
magnetic tape using an FM multiplex system. A 
maximum of 14 of these signals were demultiplexed 
and recorded on a second analog magnetic tape 
which was used for the data analysis. Due to data 
reduction hardware limitations not all desired 
combinations of signals were possible. 

Results 

Prior to applying the above analytical proce- 
dures to determine the compressor force character- 
istic from engine surge data a spectral analysis 
is performed on the available measurements. A 
typical spectral plot of compressor discharge 
pressure, P 3 J, over the entire transient being 
used is presented in Fig. 8 . Each line in the 
figure represents 2 seconds of data. The spectral 
analysis indicates that there is no rotating stall 
present during the surge phenomena. Here the 
first 16 lines (32 seconds) of the graph show the 
surge (planar) frequencies. These are a fundamen- 
tal of about 5.2 Hz with harmonics. On the 17th 
line the engine enters a non-recoverable stall 
(rotating stall); this frequency starts at 90 Hz 
and decays with engine speed. Electronic noise at 
60, 120, and 180 Hz is also seen. 

The separate occurrance of surge and rotating 
stall allows direct determination of F x , that 
is the input filters, G-j(s), in Fig. 4 are not 
required. This is quite important because some of 
the excursions in pressure and especially flow are 
very rapid. Presence of an input filter could 
soften this response to a point where the desired 
characteristic is masked, distorted, 
or becomes double valued. 


It is noted that the frequency of mode II is 
approximately two times the frequency of mode I. 


Using a rationale presented in a later sec- 
tion, output and differentiator filters with 
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corner frequencies of 200 and 400 Hz, respectively 
were chosen. 

Typical Characteristics 
Effect of Compressor Speed 

With the above filtering, compressor instan- 
taneous force characteristics, F x , were deter- 
mined at several compressor corrected speeds in 
mode I (fig. 9(a)). If it is assumed that the 
stall cell build up lag. Fig. 2, is not required 
in this period, t c = 0, then these plots are 
the compressor surge characteristics which are to 
be determined. Because of the electrical and flow 
noise on the data the results are somewhat noisy, 
however, the basic form of the characteristic is 
readily discernable. As indicated above the force 
characteristic can be converted directly to the 
compressor pressure ratio versus corrected flow 
map (fig. 9(b)). It is also of interest to com- 
pare the simple (dynamic) pressure ratio normally 
seen to these results. These plots are shown in 
Fig. 9(c) for the identical surge cycles. 

Parametric studies made with various values 
for t c during the surge instability did not 
show improvements in F m relative to F x . 

That is, accounting for the lag dynamics did not 
either reduce the "flow" noise on the character- 
istic or make the value of F during increasing 
and decreasing flow traverses coalesce better. 

If the instantaneous force characteristics 
(as presented in fig. 9(a)) are individually 
smoothed then a generalized force characteristic 
(force compressor map) may be formed, as shown in 
Fig. 10. Included on this figure for completeness 
is a characteristic determined during a mode II 
instability, the 92 percent speed line. This 
characteristic (fig. 11), (cycles 1 to 3) is some- 
what different from those of Fig. 9(a). This is 
attributed to the rapid speed change during this 
time period. 

Comparisons between adjacent mode I and mode 
II type instabilities, further into the run (when 
the speed is changing slowly), showed only small 
differences in the force characteristics. For 
N c = 88a speed, stable side data points con- 
verted from a normally measured compressor map 
using Eq. (11) are superimposed on the data of 
Fig. 10. It can be seen that the shapes of the 
stable side characteristics are the same and that 
the post-stall characteristic leaves from a region 
of high performance and force (hence pressure 
ratio) and decreases as flow approaches zero. A 
local minimum in the characteristic occurs near 
zero flow and as negative flow increases the force 
characteristic also increases. The shape of the 
characteristics is guite different from those seen 
in the literature. (*■) The post stall (positive 
flow) characteristic levels appear to be higher 
than those seen in the literature. 

A typical characteristic plot F x vs 
H:,corr shown in Figs. 12(a) and 12(b). It 
is’ annotated with letters A through G which mark 
particular events in the cycle. In Fig. 12(c), a 
set of typical variables involved in the charac- 
teristic determination are plotted against time. 
From the corresponding letter annotation it can be 
seen that the traverse from A to C (initial stall) 
and from D to F occur in a very short time 
relative to the period of one cycle of mode I. 


Typically, t/L,c is about 3.4 percent of the 
mode I period and t^,p is about 2.7 per- 
cent. These percentages will about double for a 
mode II instability. It is the rapid changes in 

W c that create a substantial W c term. This 
term then corrects the steady terms in the momen- 
tum equation (eq. (8)) to give the shape of the 
characteristic. For the same time period in the 
transient plots of P 3 T/P 2 T versus corrected 
flow and ( P 3 T/ ?2T ^ CHAR vers us corrected flow 
are presented in Figs. 12(d) and 12(e) respec- 
tively. The essential difference between Figs. 
12(d) and (e) is that (P 3 J/P 2 T) CHAR accounts 

for the W term. 

Additional Studies 

In order to examine the asumption of uniform 
flow over the length of the compressor, that is to 
determine the necessity of a spatial integration 
of the flow (eq. (2)), use of the flow at station 
3 was investigated. In the above study the flow - 
total pressure transducers were located near the 
midspan position at stations 2 and 3. A second 
analog tape was created on which hub transducers 
were used at station 2 and a forward-reverse flow 
measurement at station 3 was implied by cali- 
brating total and static pressure measurements in 
both forward and reverse flow. 

Anomalous results were found for both of 
these cases. The hub instruments appeared to 
indicate strong secondary flows and could not be 
used. Attempts to use a station 3 flow were not 
successful. The situation was further complicated 
by the fact that the total pressure was measured 
at a different circumferential location than the 
static pressure and amplitude and small phase dif- 
ferences were known to exist. Also, excessive 
instrument zero drift was seen. 

Transition to Rotating Stall 

Of particular interest in this data is the 
period of transition from the surge (mode II) 
instability to rotating stall, which occurs about 
32.5 seconds into the transient (fig. 13). Two 
methods were used to study this period of the 
transient. The first method was to use the stan- 
dard filtering (output filter corner at 200 Hz, 
differentiator filter corner at 400 Hz) and to 
examine successive sections of the characteristic 
plot. Figs. 14(a) to (c), respectively. The 
figures show the last stall followed by an 

increase in W c corr but during this traverse 
F x does not increase back to the normal charac- 
teristic levels (fig. 14(a)). Several low force 
"loops" follow (fig. 14(b) and 14(c)) with the 
final rotating stall "average" (heavily filtered) 
performance point indicated on Fig. 14(c). 

A second approach to examining the transition 
to rotating stall used progressively heavier fil- 
tering applied to the same transient. The period 
chosen began somewhat earlier (i.e., the 193rd 
cycle) and progressed into the rotating stall 
regime. The sequence shown in Figs. 15(a) to (c) 
applies successive filters each with one half the 
low pass corner frequency of the previous case. 

The force characteristic becomes distorted in 
shape and reduced in extent as the filtering is 
increased. Also, the characteristic becomes 
double valued as the filtering increases, as would 
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be expected based on the results shown in the fil- 
tering section (below). The compressor speed does 
not change significantly over the period 
considered. 

It is evident from these results that the 
"steady" rotating stall equilibrium point is sepa- 
rate and distinct from the Force (or Pressure) 
characteristics previously measured. 

Specifically, the data appears to substanti- 
ate the existance of separate surge and rotating 
stall characteristics as hypothesized by 
Moore. . 

Further study of Fig. 15(c) (and fig. 14) 
shows that the beginning of the transition starts 
from an area near the stable side of the surge 
characteristic (near a minimum P 37 and maximum 

W q C o rr ) and progresses to the left (reduced 
flow) and settles on the rotating stall charac- 
teristic. This is also seen from the variables 
plotted against time in Fig. 13, note beginning of 
rotating stall on P2T- Thus, it is conjectured 
that the combination surge-rotating stall charac- 
teristic has a form similar to that shown in Fig. 
16. This form assumes some sort of a bifurcation 
point connecting the two characteristics. If this 
conjecture can be substantiated with additional 
data, the the important research question then 
becomes one of determining the physical (and hence 
mathematical) criterion which causes switching 
from one characteristic to the other. It is 
believed that this criterion will involve a 
relationship between the internal forces com- 
posing F x . It is for this reason that the 
characteristics in this paper are presented as 
forces rather than pressure ratios. How the 
Grietzer B criterion^) fits into this picture 
is not clear but it is suspected that other issues 
which have not before been considered may be 
important. One such issue is the amplitude of the 
surge cycle, since if a bifurcation point does 
indeed exist switching from the surge characteris- 
tic to the rotating stall characteristic can only 
take place if the amplitude of the surge oscil- 
lation (on F x versus W c corr map) reaches the 
bifurcation point . 

Indeed, if this truly three-dimensional prob- 
lem of compressor stability is to be grossly sim- 
plified into a two oscillator model, one for surge 
and one for rotating stall, two separate and dis- 
tinct criteria may be required. That is, to 
explain those cases where rotating stall oscil- 
lations occur stably over part of the period of 
the surge cycles would require an on-off criteria 
for the rotating stall oscillator. While to cover 
the transition from surge to rotating stall (or 
reverse) would require a second criterion. Such 
conclusions, of course, can not be drawn on the 
basis of a single study and require consideration 
of dynamic data from many compressors. 

Combining the rotating stall point (at 
N c = 84 percent) with the generalized surge 
force characteristic (at N c = 84 percent) yields 
the final combined characteristic map (fig. 17). 


Filtering 

Filtering is a crucial step in the determina- 
tion of compressor characteristics from surge 
data. As mentioned above, filters are required to 
remove; flow and electrical noise, effects of 
rotating stall if present, and noise from the dif- 
ferentiation process. Because of the rapidity of 
certain parts of the surge transient (high fre- 
quency content) only very light filtering may be 
used. The particularly sensitive parts of the 
transient are rapid movements indicated on Fig. 
12(b) from A to C and from D to F. Referring to 
F x vs time (Fig. 12(c)), these segments appear 
as very short term transients, and are difficult 
to separate from background noise. 

The basic phi losophy. f or selecting filter 
time constants (hence corner frequencies) was as 
follows. The filter was used with the lowest 
corner frequency which met the following criteria; 
(a) further increasing the corner frequency did 
not affect the A to B and D to F traverses in Fig. 
12 (b) materially, and (b) further decrease of the 
corner frequency would affect the basic shape of 
the characteristic or open the gap between the A 
to B and D to E traverses. This process was 
applied to the differentiator filter (for F x 
calculation) first, then to the output filters. 
Keeping the output filters in matched pairs 
appears to allow that corner frequency to be less 
than that of the differentiator filter. This 
resulted in corner frequencies of 400 and 200 Hz 
for the differentiator and output filters, 
respectively. The filter required for the dif- 
ferentiator used to calculate F m , was set by 
allowing t c = 0 , and increasing the filter 

corner frequency until F m vs W c>corr 

approximated F x vs W c corr . Then t c was 
set to its desired (research) values. This 
yielded a corner frequency of 2400 Hz. 

The high corner frequencies chosen allowed a 
significant level of noise on the plotted results, 
however, this was felt to be justified to cor- 
rectly identify underlying characteristic. The 
effects of moving the output and differentiator 
filter corner frequencies relative to the selected 
values is shown in Fig. 18. Reducing the corner 
frequencies tends to spread the rapid traverses 
indicated above making the characteristic appear 
to be double-valued. Increasing the corner fre- 
quency increases the noise without improving the 
definition of the underlying characteristic. 

Summary and Conclusions 

A technique for extracting the in-stall 
pumping characteristics for an axial flow compres- 
sor from engine surge data has been developed. 

The basic concept appears to be workable and is 
premised on uniform flow (time varying) over the 
compressor. The technique has been applied to a 
multi-stage compressor in a turbofan engine. 

The manner in which the data are filtered has 
been found to be critical in the successful appli- 
cation of the method. The separate occurrance of 
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APPENDIX A 


surge and rotating stall in the data analyzed sim- 
plified the filtering problem, in as much as input 
filters were not required. Careful selection of 
the filters required for the approximate differ- 
entiations was also necessary. Excessive fil- 
tering was found to produce distorted and double- 
valued surge characteristics. Insufficient 
filtering left excessive noise on the resulting 
characteristic. 

Two instability modes were found in the 
data. These were a classical surge mode and a 
second instability at about half the surge fre- 
quency thought to be associated with the engine 
control. Characteristics determined from both 
modes, when compressor speed was near constant, 
were found to be indistinguishable. 

Consideration of the transition from surge 
into rotating stall indicates the existence of 
separate surge and rotating stall characteris- 
tics. It appears that these characteristics meet 
at a bifurcation point located on the stable side 
of the surge characteristic. 

The stall cell build-up lag dynamics were 
considered and were found not to improve (elimi- 
nate dynamics from) the surge characteristic. The 
lag thus, may not to apply during the surge pro- 
cess. Accounting for the build-up lag did not 
yield improved results in the transition from 
surge to rotating stall. 


Expressions are developed to allow further 
investigations into the nature of the components 
of the net compressor force. Use of these 
requires an estimate of the compressor frictional 
behavior. 


Finally, the characteristics determined by 
the procedure are felt to be representative of the 
true characteristics. However, further refine- 
ments in the areas of instrumentation type, place- 
ment, quantity, noise, and drift levels an 
improved data acquisition and filtering technique 
should allow significant reduction of noise and 
improved definition of the characteristics and 
their interrelationship. 


ALTERNATE MOMENTUM EQUATION 

Rewriting Eq. (5) as: 

2(P t - P S )A = P AV2 (Al) 

and substituting Eq. (Al) into Eq. (4) and solving 
for F x yields: 

F y = (2 P^y - P^j) A-j - (2 P 7T - P, c ) A, 


2T ” 2S' 


as an alternative for Eq. (8). 
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Fig. 2. - One dimensional surge model. 
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Fig. 3. - Simplified block diagram for determination of F x and F m . 









Fig. 4. - Block diagram for practical determination of forces F x and F m . 



Fig. 5. - Instantaneous force components. 













3S” 3S, nom 



TIME 


Fig. 6. - Non recoverable stall - transient. 
Fuel pulse at 92% N c 
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Fig. 7. - Interstage static pressures. 
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Fig. 8. - Power spectral analysis versus time of P^j 
surge into rotating stall at 2 sec/line. 


1.2 



(a -3) N c ■ 84. 7%. 


(a) Compressor force characteristic. 

Fig. 9. - Effect of compressor speed. F 0 = 200 Hz*. F^ = 400 Hz. 
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(b-3) N c = 84.7%. 

or pressure ratio characteristics. 
Fig. 9. - Continued. 




w 2 c, corr 
<c-3) N c = 84. 7%. 


(c) Compressor pressure ratios. 
Fig. 9. - Concluded. 
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Fig. 10. - Compressor generalized force characteristic. F o = 200 Hz; 
F d = 400 Hz. 



Fig. 11. - Force characteristic in mode II instability. N c = 92%; 
F 0 = 200 Hz; F d = 400 Hz. 





(a) Force characteristic. 


Fig. 12. - Typical transient. N c = 86. 4%; F 0 = 200 Hz; F d = 400 Hz; 
cycle 22. 



(b) Annotated force characteristic. 
Fig. 12. - Continued. 
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(d) Pressure ratio. 


Fig. 12. - Continued. 



(e) Pressure ratio characteristic. 
Fig. 12. - Concluded. 
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(b> Time segment 2. 
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(c) Time segment 3. 

Fig. 14 - Transition to rotating stall - progressive time segments. 
N c = 84%; F 0 = 200 Hz; F d = 400 Hz. 












<c) F 0 = 25 Hz; F d = 50 Hz. 

Fig. 15. - Transition to rotating stall - progressive filtering, 
N c = 84%. 





Fig. 16. - Combined surge - rotating stall characteristic. 
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Fig. 17. - Compressor generalized force characteristic. F 0 = 200 Hz ; 
F d = 400 Hz. 



(d) Decreasing differentiator filter corner frequency, F 0 *= (e) Decreasing output filter corner frequencies; F 0 = 50 Hz; 

200 Hz; F d = 50 Hz. F d = 400 Hz. 


Fig. 18. - Filter study. N c = 86.6%; cycle 14. 
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